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ABSTRACT 

Consider  an  n  component  reliability  system  having  the 
property  that  at  any  time  each  of  its  components  is  either 
up  (i.e.,  working)  or  down  (i.e.,  being  repaired).  Each 
component  acts  independently  and  we  suppose  that  each  time 
the  ith  component  goes  up  it  remains  up  for  an  exoonentially 
distributed  time  having  mean  ,  and  each  time  it  goes 
down  it  remains  down  for  an  exponentially  distributed  time 
having  mean  v^  .  We  further  suppose  that  whether  or  not 
the  system  itself  is  up  at  any  time  depends  only  on  which 
components  are  up  at  that  time.  We  are  interested  in  the 
distribution  of  the  time  of  first  system  failure  when  all 
components  are  initially  up  at  time  zero.  In  Section  2  we 
show  that  this  distribution  has  the  NBU  (i.e.,  new  better 
than  used)  property,  and  in  Section  3  we  make  use  of  this 
and  other  results  to  obtain  a  lower  bound  to  the  mean  time 
until  first  system  failure. 
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1.  INTRODUCTION  AND  SUMMARY 

Consider  an  n  component  reliability  system  having  the  property  that  at  any 
time  each  of  its  components  is  either  up  (i.e.,  working)  or  down  (i.e.,  being 
repaired).  Each  component  acts  independently  and  we  suppose  that  each  time  the 
ith  component  goes  up  it  remains  up  for  an  exponentially  distributed  time  having 
mean  ,  and  each  time  it  goes  down  it  remains  down  for  an  exponentially  distri¬ 
buted  time  having  mean  v^  . 

We  suppose  that  whether  or  not  the  system  itself  is  up  at  any  time  depends 
only  on  which  components  are  up  at  that  time.  In  particular,  letting 


1  if  the  ith  component  is  up  at  time  t 


Xi(t)  = 

(0  otherwise 


f  ‘ 

(0  o 


if  the  system  is  up  at  time  t 


therwise 


we  suppose  that  there  exists  a  nondecreasing  binary  function  d>  such  that 


X(t)  =  4> (X.  (t) ,  ...,  X  (t))  . 
1  n 


The  function  <}>  is  called  the  structure  function  o’'  the  system. 

The  above  system,  with  the  added  generality  of  arbitrary,  rather  than 
exponential,  component  uptime  and  downtime  distributions  was  considered  in  Ross  [5] 
where  such  quantities  of  interest  as  (i)  the  average  system  failure  rate;  (ii)  the 


average  length  of  system  uptime;  and  (iii)  the  average  length  of  system  downtime, 
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were  computed.  In  the  present  paper,  we  are  interested  in  the  distribution  of  the 
time  of  first  system  failure  when  all  components  are  initially  up  at  time  zero.  In 
Section  2  we  show  that  this  distribution  has  the  NBU  (i.e.,  new  better  than  used) 
property,  and  in  Section  3  we  make  use  of  this  and  other  results  to  obtain  a  lower 
bound  to  the  mean  time  until  first  system  failure. 


SJWWT"’ 
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2.  SOME  PRELIMINARIES  CONCERNING  NBU  DISTRIBUTIONS 


We  say  that  a  lifetime  distribution  F  is  NBU  (new  better  than  used)  if 


F(t  +  s) 


<_  F(s)  for  all  s  ,  t  >_  C 


In  words,  F  is  N1U  if  the  probability  that  a  t  year  old  item  survives  an 
additional  s  years  is  less  than  or  equal  to  the  probability  that  a  new  item 
survives  s  years,  for  all  s  ,  t  . 


Proposition  1; 

For  an  arbitrary  structure  .  Assuming  that  all  components  are  initially 
up,  the  time  until  the  first  system  failure  has  an  NBU  distribution. 

Proofs 

We  shall  prove  the  proposition  by  first  showing,  by  induction,  that  the  time 
to  first  system  failure  will  have  a  stochastically  larger  distribution  if  all 
components  are  initially  up  as  opposed  to  n  -  k  of  them  being  initially  up  and 
the  other  k  initially  down.  To  prove  it  for  k  =  1  ,  suppose  that  components  2 
through  n  are  initia  l  ly  up  and  suppose  that  there  are  two  identical  components 
that  can  be  used  as  component  1  — namely  1'  which  is  initially  up  or  1'  '  which 
is  initially  down.  To  show  that  the  system  will  have  a  stochastically  larger  time 
to  first  failure  if  we  use  1'  as  opposed  to  1''  ,  let  us  condition  on  the  first 
time  t  such  that 


4>(1,X2(t),  ....  Xn(x))  =  1  ,  <K0,X2(t) . VT))  =  0  * 


Now,  given  that  t  =  t  ,  the  probability  that  the  first  system  failure  will  occur 
before  (or  at)  time  x  is  the  same  irregardless  of  whether  1'  or  1''  is  used 
when  x  <  t  .  On  the  other  hand,  if  t  =  t  and  the  system  has  not  failed  prior  to 


time  t  ,  then  the  probability  that  it  fails  at  time  t  is  equal  to  the  probability 
that  component  1  is  down  at  t  .  However,  given  that  t  ■  t  ,  it  follows  from 
the  theory  of  two-state  Markov  processes  that  the  probability  that  the  first 
component  is  up  at  time  t  is  greater  if  it  is  initially  up  as  opposed  to  being 

/  vi 


initially  down  (the  probabilities  being 


h  +  vi  "i  +  vi 


-  (l* .  +V .  )  t 

e  and 


■<ui+v1> 


"i  +  vi 


“1 +  vi 


') 


Therefore,  given  i  =  t  ,  the  system  is  more  likely 


to  experience  its  first  failure  by  time  t  if  l'1  is  used.  Furthermore,  by  the 
lack  of  memory  of  the  exponential  distribution,  it  follows  that  if  the  system  has 
not  failed  by  time  t  (and  thus  component  1  is  up  at  time  t)  then  the  distribu¬ 
tion  of  the  additional  time  until  system  failure  is  the  same  irregardless  of  which 
component — l1  or  1  *  * — was  used  as  component  1  .  Hence,  given  t  ,  the  system 
will  have  a  stochastically  larger  time  to  first  failure  if  1* — as  opposed  to  1M — 
is  used.  Taking  the  expectation  of  x  (i.e.,  uncondicioning) ,  then  yields  the 
result  when  k  =  1  . 

The  general  case  is  now  easily  established  for  if  k  of  the  components  are 
initially  down  then  by  concentrating  on  one  of  these  components  the  same  argument 
as  the  one  presented  for  the  case  k  =  1  shows  that  it  would  have  been  better  (in 
a  stochastic  sense)  if  this  component  was  initially  up.  But  by  the  induction 
hypothesis  all  components  initially  up  is  better  than  having  any  set  of  k  -  1 
initially  down  and  thus  the  induction  is  complete. 

The  proposition  now  follows  for  if  the  system  is  up  at  time  t  we  see,  by 


conditioning  on  the  number  of  components  up  at  time  t  ,  that  the  system  has  less 
chance  of  surviving  an  additional  s  years  than  does  a  system  that  has  all  of  its 
components  up  (i.e.,  a  new  system) .  |  | 


f 
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Remark ; 

If  all  components  are  not  initially  up,  then  Proposition  1  is  no  longer  valid. 
For  instance,  consider  a  parallel  system  of  two  components  (i.e.,  the  system  is 
down  when  both  components  are  down)  which  have  identical  exponential  uptime  and 
downtime  distributions,  and  suppose  that  one  of  the  components  is  up  at  time  zero 
and  the  other  is  down.  Then  if  the  system  has  not  yet  failed  by  time  t  it  follows 
that  the  remaining  time  to  failure  is  a  mixture  of  two  distributions  one  of  which 
is  the  distribution  of  the  time  from  zero  to  the  first  failure  while  the  other  one 
is  (strictly)  stochastically  larger.  (The  first  distribution  comes  into  effect  if 
only  one  of  the  components  is  up  at  time  t  and  the  other  if  both  are  up  at  time 
t  .)  Hence,  the  distribution  of  additional  time  to  failure  would  be  stochastically 
larger  than  the  distribution  of  time  to  first  failure.  In  fact,  it  follows  from 
results  presented  by  Keilson  [3]  that  the  distribution  of  time  to  first  failure  is 
in  this  example,  a  mixture  of  exponential  distributions,  which  not  only  is  not  an 
NBU  distribution  but  is  a  DFR  (decreasing  failure  rate)  distribution  (see  [5]  for 
definition. j  | 

A  second  result  about  NBU  distributions  that  we  shall  use  is  embodied  in  the 
following  proposition. 

Proposition  2: 

Let  T, ,T- ,  ....  T  be  nonnegative  independent  random  variables,  each  having 
i  i.  m 

an  NBU  distribution,  and  let  =  E [ T^ ]  .  Then 

E[Min  (T1#  ....  Tm)]  >  |j  1/p.j  . 

Proof ! 

Let  X  be  an  NBU  distribution  having  distribution  F  and  mean  and  let 
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Y  be  an  arbitrary  nonnegative  random  variable  having  distribution  G  and  mean 
Wy  ,  and  assume  that  X  and  Y  are  independent.  Now, 


Also, 


CO 

E[Min  (X,Y) ]  -  J  E[Min  (X,t)]dG(t)  . 


E[Min  (X, t)  ]  =  xdF(x)  +  t(l  -  F(t)) 


W 


*  V  -  I  xdF(x)  +  t(l  -  F(t) ) 


Ux  -  E[X  |  X  >  t ]  (1  -  F(t) )  +  t(l  -  F(t) ) 


>  Px  ~  (t  +  Wx)(l  "  F(t) )  +  t(l  -  F(t) ) 


F(r)u 


when  we  have  used  the  fact  that  F  NBU  implies  that  F. [X  |  X  >  t]  <_  t  +  y^ 
Hence, 


E[Min  (X,Y) ]  >  ^  J  F(t)dG(t) 


=  y  P{X  <  Y} 
x  — 


Now  let  T  ,  i  -  1,  ....  m  be  independent  NBU  random  variables  with  means 
y  ^  ,  i  -  1,  . ..,  ra  .  From  Equation  (1)  we  obtain,  for  each  i  , 


:[Min  (1^,  ....  T  )]  =  E[Min  (Ti  ,  Min  (T^  ...,  ...,  Tm))] 


>  y,P{T,  =  Min  (T  ,  ....  T  ) } 
—  1  i  i  m 
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Hence,  dividing  the  above  by  and  then  summing  Dver  i  yields 


(jt  (Tj . V1  PtTi 


Min  (T.  ,  ....  T  )} 
i  m 


and  the  proof  is  complete. 


Remark : 

Proposition  2  states  that  the  mean  life  of  a  series  system  of  independent 
NBU  components  is  greater  than  or  equal  to  the  mean  life  of  a  series  system  of 
independent  exponential  components  having  the  same  set  of  means.  (A  reverse 
inequality  can  be  proven,  in  the  same  way,  for  parallel  systems.)  This  result  had 
previously  been  proven  (see  [4])  under  the  stronger  (than  NBU)  assumption  that  each 
component  lifetime  is  IFRA,  where  we  say  that  the  lifetime  distribution  F  is 
IFRA  (increasing  failure  rate  on  the  average)  if  lo8  is  nondecreasing 

in  t  .  Further  we  see  from  the  proof  of  this  proposition  that  the  condition  that 
the  are  NBU  can  be  reduced  to  the  condition  that  E[T^  |  >  t]  <_  t  +  E[T^] 

for  all  t  .  This  weaker  than  NBU  assumption  is  referred  to  in  the  literature  as 
assuming  that  T^  is  NBUE  (new  better  than  used  in  expectation). 


After  formulating  and  proving  Proposition  2,  the  author  discovered  that  the  result 
is  contained  as  an  exercise  in  a  forthcoming  text  on  reliability  theory  by 
R.  Barlow  and  F.  Proschan.  However,  f.neir  indicated  method  of  proof  is  different 
from  the  one  employed  in  the  present  paper. 
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3.  MEAN  TIME  TO  FIRST  SYSTEM  FAILURE 

One  way  of  characterizing  the  structure  function  $  of  the  reliability  system 
is  in  terms  of  its  minimal  cut  sets.  We  say  that  the  subset  C  of  the  set  of 
components  is  a  cut  set  if  the  system  is  necessarily  down  when  all  components  in 
this  subset  are  down.  In  addition,  we  say  that  it  is  a  minimal  cut  set  if  it  does 
not  (properly)  contain  any  other  cut  sets.  Let  us  denote  by  C^,  ...»  Cg  ,  the 
minimal  cut  sets  of  the  structure  $  .  Since 

4>(x^,  ....  xn)  -  Min  (o^,  ....  ag) 


where 


!0  if  all  components  of  are  down 

1  otherwise 

it  follows  that  knowing  the  mirimal  cut  sets  is  equivalent  to  knowing  the 
structure  <J>  . 

If,  in  our  reliability  system,  we  let  denote  the  first  time  that  all 

components  in  the  ith  cut  set  are  down,  then  T  ,  the  time  to  first  system  failure, 
is  given  by 

T  =  Min  (T, ,  ...,  T  )  . 

As  the  so  defined  will  not,  in  general,  be  independent  (they  will  be 

independent  only  if  no  two  minimal  cut  sets  overlap)  we  cannot  immediately  apply 
Propositions  1  and  2  to  obtain  a  lower  bound  for  E(T)  in  terms  of  the  E[T^]  . 
However,  it  has  been  shown  by  Esary  and  Proschan  [2],  using  the  concept  of  associa¬ 
tion  of  random  variables,  that  if  the  component  uptime  and  downtime  distributions 


are  exponential  then 
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8 

p{t  >  t}  >_  n  p{t  >  t)  . 

i-l  1 


Using  this  result,  now  yields 


Proposition  3: 

Assuming  that  each  component  is  Initially  up,  then 


E(T]  > 


1/EITJ 


Proof : 


E[T] 


P{T  >  t }dt 


IlP{Tj  >  t}dt  by  Esary  and  Proschan  [ 2 ] 


> 


l/E(Ti> 


by  Proposition  2. | 


Thus,  it  remains  to  determine  ECT^  .  Now  is  the  first  time  that  all 

components  of  the  ith  minimal  cut  set  are  down,  and  so  considering  only  those 

components  in  we  see  that  determining  E(T^)  is  equivalent  to  determining 

the  mean  life  of  a  parallel  system  when  all  components  are  initially  up.  So  let 
us  consider  the  parallel  system  of  the  (say)  r  components  of  the  minimal  cut  set 
C  and  suppose  that  these  r  components  have  mean  exponential  uptimes  p^,  ...,  pf 
and  mean  exponential  downtimes  v^,  ...,  .  Now  it  has  been  shown  by  Brown  [1] 

that  the  mean  time  until  all  components  of  C  are  down  is  equal  to 
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k-!  l1<i2<...<ikM-1  ij 


k  Ui 


Thus,  we  can  use  Brown's  formula  above  to  compute  E[T^]  for  each  minimal  cut  set 
and  then  use  Proposition  3  to  obtain  a  lower  bound  on  the  mean  time  to  first  system 
failure.  These  computations  can  be  easily  effected  by  use  of  a  computer. 
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